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Abstract 

In this paper we give a first attempt to define and study stable dis- 
tributions with respect to the weak generahzed convolution, focusing 
our attention on the symmetric weakly stable distribution. As in the 
case of the classical convolution, characterization of distributions sta- 
ble in the sense of the weak generalized convolution depends on solving 
some functional equations in the class of characteristic functions. 
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1 Introduction 



The investigations of weakly stable random variables started in the seven- 
ties in the papers of Kucharczak and Urbanik (see [TT], El])- Later a series 
of papers on weakly stable distributions written by Urbanik, Kucharczak 
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and Vol'kovich was published (see e.g. [121 EH EOl EI]). Recently a paper 
written by Misiewicz, Oleszkiewicz and Urbanik (see [I7]) appeared, where 
one can find a full characterization of weakly stable distributions with non- 
trivial discrete part, and a substantial attempt to characterize weakly stable 
distributions in the general case. 

In stochastic modeling of real processes, using independent random vari- 
ables or Gaussian processes in a variety of constructions turned out to be not 
sufficient or adequate. Multidimensional stable distributions have nice lin- 
ear properties and enable more complicated structures of dependencies, thus 
recently the role of stable processes in stochastic modeling is growing. On 
the other hand stable distributions are very difficult in calculations. However 
there are some efficient techniques for their computer simulation. 

In this situation weakly stable distributions and processes seem to be 
good candidates for use in stochastic modeling. They extend Feller's idea of 
subordinated processes. They have nice linear properties, namely: if (Xj) 
is a sequence of independent identically distributed random vectors with a 
weakly stable distribution, then every linear combination J2 cn^i has the 
same distribution as Xi ■ 6* for some random variable 9 independent of Xi. 
This condition holds not only when (oj) is a sequence of real numbers, but 
also when (a,) is a sequence of random variables such that (oj) and (Xj) 
are independent. This means that dependence structure of the linear com- 
bination ^ OiXj and dependence structure of the random vector Xi are the 
same, and the sequence (oj) is responsible only for the radial behavior. More- 
over, weak stability is preserved under taking linear operators, projections 
or functionals. On the other hand, radial properties of distribution can be 
arbitrarily defined by choosing a proper random variable 6 independent of 
Xi and considering the distribution of 6* ■ Xi. Similar properties of tempered 
stable distributions (see e.g. |20]) are the reason why they are so important 
now in statistical physics modelling turbulence, or in mathematical finance 
for modelling stochastic volatility. 

In this paper we develop the idea of distributions stable with respect to a 
generalized convolution defined by weakly stable variables. In this construc- 
tion the weakly stable variable plays the role of a catalyst in the presence of 
which the underlining process can develop. This weakly stable variable can 
be also treated as a filter, so that we observe the original process only by its 
filtered values. In this sense, we want to characterize distributions which are 
stable after filtering. 

By V(K) we denote the set of all probability measures on a separable 
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Banach space E (with dual E*). For simplicity we write V for the set of 
all probability measures on M. Moreover, V+ denotes the set of probability 
measures on [0, cxd). The symbol $ stands for the set of all characteristic 
functions on M, whereas (5x denotes the probability measure concentrated at 
the point x G E. If a sequence of probability measures (An)neN converges 
weakly to a probability measure A, we write A„ ^ A. 

For a random vector (or random variable) X, we write vC(X) for the 
distribution of X. For A e A = C{e), we write |A| = C{\e\). If /i = 
£(X) G 'P(E), then the characteristic function /i: E* — C of the measure ix 
(of the random vector X) is defined by 



For random vectors X, Y, we write X = Y for >C(X) = C{Y). If X and Y 
are independent random vectors then £(X + Y) is the convolution of >C(X) 
and C{Y), denoted by C{X) * C(Y). 

For t G M, a rescaling operator T^: V{K) — > 'P(E) is defined as follows: 



It is easy to see that if /i = C(X.) then T^/i = £(tX). The scale mixture fio X 
of the measure /i G V(E) with respect to the measure \ eV is defined by 



If fi = C(X.) and A = C{9) with X and 9 independent, then fio X = C(K9). 

A random vector X with the distribution on a real separable Banach 
space E is weakly stable iff 



where X' is an independent copy of X and the random variable 9 is inde- 
pendent of X. It was proved in [I7] that X is weakly stable if and only 




exp {i < C,,x >} ii{dx). 





Va,6GM3e aX + 6X' = X^, 



{A) 



if 



\/ 9^,92 ^9 9iX + 92X' = X9, 
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where 9i, 62 are real random variables such that 61,62, X, X' are independent 
and the random variable 9 is independent of X. In the language of probability 
measures, the condition (B) can be written as 

\/Xi,X2eV BXeV (/i o Ai) * (/i o Ai) = o A, (C) 

for A, Aj the distributions of 9,9i, i = 1,2. It was shown in [17] that the 
measure A is uniquely determined if the measure /x is not symmetric. For a 
symmetric measure /i, we only have uniqueness of the measure |A|. 

The best known examples of weakly stable random vectors are symmet- 
ric stable vectors, and in this case the random variable 9 appearing in the 
condition (A) is a constant, 9 = {\a\°' + for some a G (0,2]. Another 

family of weakly stable distributions consists of uniform distributions Un on 
unit spheres Sn-i C M" and their lower-dimensional projections. 

It was shown in \17\ that if a weakly stable distribution /i contains a 
discrete part, then it is discrete and either /i = (5o, or /x = ^6a + |5-a for 
some a G E \ {0}. From now on we will assume that the considered weakly 
stable measure /i is non-trivial in the sense that it is not discrete. 

We can now define a generalized weak convolution © = for any non- 
trivial weakly stable measure /x. 

Definition 1 Let X be a non-trivial random vector with the weakly stable 
distribution fi. The weak generalized convolution © of measures Ai, X2 E V 
is defined by 

^ ^ ^ _ i X if fi is not symmetric, 
^ ^ \ l-^l if l^ is symmetric, 

where X G V is such that (/i o Ai) * (n o X2) = jj, o X. For two independent 
random variables 9i and 6*2 with distributions Ai and A2 respectively, the weak 
generalized sum 9i © 92 is the random variable defined by 

9iX + 92X' = X (01 © 92) , 

where 9i, 92, X, X' are independent, the random variable 0i©02 is independent 
ofl^ andC{9i®92) = Ai © A2. 

The operation © in "P is commutative and associative. Moreover, as 
shown in [18], the following conditions hold: 

(i) the measure So is the unit element, i.e. © A = A for all A G P if /i is not 
symmetric and (5o © A = |A| for all A G P if /i is symmetric; 
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(ii) {p\i + q\2)®\ = p(Ai©A)+g(A2©A), whenever A, Ai.As G P andp ^ 0, 
q^O,p + q = l (linearity); 

(iii) (TaAi)©(TaA2) = Ta(Ai©A2) for any Ai, A2 G P and a > (homogeneity); 

(iv) if A„ —>■ Ao then A„ © A ^ Aq © A for all A G P (continuity). 

The idea of generalized convolutions has been extensively studied after it 
was introduced by K. Urbanik in 1964 [25]. The definition proposed by K. 
Urbanik is as follows: 

A commutative and associative binary operation o : V+ x V+ 
is called a generalized convolution if it satisfies conditions (i)-i-(iv) with © 
replaced by o and the following condition holds: 

(v) there exists a sequence (c„)„gN of positive numbers such that the sequence 
{Tc^5\^)nm weakly converges to a measure different from 5q. 

The first, but not the most important difference between the definition 
of generalized convolution o given by K. Urbanik and the definition of weak 
generalized convolution © is the domain, i.e. ® : V ^ V. This implies, in 
particular, that most of the methods used in studying generalized convolu- 
tions cannot be directly applied for weak generalized convolutions. Another 
difference is that the weak generalized convolution have properties (i), (ii), 
(iii) and (iv), but, in general, not (v). In spite of this disadvantage, we do not 
have to assume that the algebra {V, ©^) is regular, i.e. that there exists a 
non-trivial homomorphism of (P, ©^) into a complex field. The assumption 
that there exists a non-trivial homomorphism from {V+, o) into the positive 
half-line was crucial in studying generalized convolutions. In the case of the 
weak generalized convolution ©^ we have that for every ,^ G E* the formula 



defines a homomorphism of (P, ©^) into a complex field. Moreover, if the 
weakly stable measure is non-trivial, then there exists ^ G E* such that /i^ 
is non-trivial. 

Talking about the homomorphism and treating the set V as an algebra 
we underline here that V is equipped with the generalized convolution ©^ as 
a binary operation and with rescaling measures operator Tf, t G M, which can 
be treated as multiplication by scalars. Moreover, we see that convex linear 
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combinations of probability measures are also probability measures, thus V 
can be treated subset of some linear space. 

In this paper we consider only symmetric weakly stable distributions /i, 
thus the functions /i^, ^ G E*, and /i are real- valued. Moreover, for symmetric 
we know that //oA = /i.o|A|,sowe can restrict our attention to the set 
{V+,®^) instead of {V,®^). 

The paper is organized as follows. In section 2 solve a functional equa- 
tion, which will be needed later. In section 3 we describe the problem of 
characterizing stable distributions with respect to the generalized convolu- 
tion The description of strictly stable distributions in this sense is given 
in subsection 3.1. In the following subsections we discuss the possibility of 
having occasionally some strange behavior of the function d. The last subsec- 
tion contains discussion of the general case of stable, but not strictly stable, 
distributions in the sense of the generalized convolution 

2 Functional equation 

Let be the set of all continuous functions /: [0, oo) — > R such that /(O) = 
and is an isolated point of the set {x e M: f{x) — 0}. 

Theorem 1 Let a,b E (0, oo) and let f & T be a solution of the functional 
equation 

f{t)^f{at) + f{bt). (1) 

Then a,b E (0, 1) and there exist p > and a continuous function H: (0, oo) — > 
M \ {0} such that H{t) = H{at) = H{bt) for every t > 0, and 

f{t)=fH{t), t>0. 

Moreover, if In a /In b is irrational then the function H is constant. 

The proof of Theorem 1 is based on a series of lemmas. The first one can 
be easily proved by mathematical induction. 
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Lemma 1 Let a,b e (0, oo) and let f: [0, oo) — > M 6e a solution of equation 
(1). Then for every n E N and every t eM vje have 



k=0 

□ 



Lemma 2 Let a,b E (0, oo). // (1) has a solution in the class T then 
a,be (0,1). 

Proof. Let / G be a solution of (1). Without loss of generality we may 
assume that b and /(x) > for every x e (0,m) with some m > 0. Then 
for every t G (0,m) we have ba~^t G (0,m) and, consequently, 

fit) = f{aa~H) = f{a-h) - f{ba~-H) < fia^H). 

If a > 1 this and mathematical induction would imply that f{t) ^ f{a~^t) 
for every n G N, and by the continuity of / we would have /(t) = for every 
t G (0,m). If a = 1 then equation (1) imply that f{bt) = for each t > 0. In 
both cases f T which contradicts our assumption. Therefore a < 1. □ 

Lemma 3 Let a,b E {0,oo) and let f E be a solution of equation (1). 
Then f(t) = if and only ift = 0. 

Proof. Assume, for instance, that for some m > we have /(x) > for 
every x E (0, u). By Lemma 2 we know that a,b E (0, 1), thus for every fixed 
t > we can find n G N large enough to have 

a^b'^-h E {0,u), A; = 0, l,...,n. 

Then, by Lemma 1, we have 

which ends the proof. □ 

The proofs of Lemmas 4-6 below take pattern of some ideas from |^ (see 
also tU). 
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Lemma 4 Let i,j e Z. For every a > and (3 e (0,2q;) there exists 
no ^ max{|i|, such that for every n ^ 2no 



k+i ( k-i J^^U 



A; e {no, ...,n- no}. 



Proof. Let a > and /? e (0, 2q;) be fixed. Assume first that i + j ^ 0. 
Then at least one of the numbers i.j, for example i, is nonnegative. There 
exists no e N such that no ^ max{|i|, |j|} and 



1 - 



2a 



no + 1 / V n-o + i. 
Fix integers n and k such that n ^ 2no and no ^ A; ^ n — no- Then we have 
n — i — f 



A; + i / V k — i 



Wn + 



a (n + 1) . . . (n + i + j) 



1 



k+1 



\J\ 



n 



where 



k + i J \ n — k + 1 

(n + 1) . . . (n + i + j) ■ (n - A;)! 
(A; + l)...(A; + i) • (n- A; + j)!' 



j + 1) . . . n 

n-k+ \j\ 



This implies that 

^n + i + j 



A; 



+ a' 



n 



k — i 



n 



Wn + 



Wr. 



1 - 



no + 1 



1 - 



\J\ 



no + 1 



4 Wr, 



n 



2 4/3, 



n 



/5, 



where the last inequality follows from the fact that the function (0, oo) 3 x ^ 
X + /?^/4,T attains its minimal value at a; = (3/2. In the case i + j < it is 
enougii to replace i,j,a,l3 by —i, —j, 1/q;,/3/q;^, respectively, in the previous 
reasoning. □ 
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Lemma 5 Let a e (0, oo), h e (0, 1) and let f: [0,oo) R be a solution 
of (1) continuous at zero and such that /(O) = 0. Then for every k & N and 
every t > we have 



n=k ^ ^ 



In particular, for every k &N, and every t > 



lim f {a''b''-H) =0. 



Proof. Let t > 0. By mathematical induction we have that for every m e N 



m—l 



fit) = J2 /(«&"^) + /(&™^) ^ 



n=0 



n=0 



since b e (0, 1) and / is continuous at zero, so fiV^t) — > for m — > cxo. This 
means that we have proved the required equality in the case A; = 0. This 
equality will be used in the next step of the proof. 

Now let /c e N be fixed and assume that for each i > 



n=k-l ^ ^ 



Then we have 

oo 



n=k-l ^ ^ n=k ^ ^ 

n=k ^ ^ 3=0 j=0 n=k ^ ^ 



CXJ CXJ / ^ \ CXJ 

EE = E 



E 

,n=k 



n — 1 

k - 1 



which ends the proof. 



□ 
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Lemma 6 Let a,b e (0, oo) and let f e be a solution of (1). Then for 
every i,jEZ 

f{tf ^ f (a-^b-H) f (a'l^t) , t>0. 

Proof. Lemma 2 yields a, 6 G (0, 1). By Lemma 3 we may assume, without 
loss of generality, that / is positive. Let i,j G Z and a > 0, /5 G (0,2a) be 
fixed. Choose no G N for which the assertion of Lemma 4 holds. Let t > 
and e > 0. By Lemma 5 there exists n ^ 2no such that 



n 



for every /c G {0, 1, . . . , no — 1}. Then 

Q / ^ A: G {0, . . . , no - 1} U {n - no + 1, . . . , n}. 

Now, by Lemma 1, we obtain 

/ {a-%~H) + a^f {a}Vt) 

n+i+j , . .s n—i—j 



'-t-'-t- J / _|_ • _|_ -X / _ ■ _ A 

k=0 ^ ^ fc=0 ^ 

m— rid ^ ' rr)— n ^ ' 



m=—i 
n— no 

:^ , 

m=n{) ^ m=0 

- e. 



Since /3 was an arbitrary number from the interval (0, 2a) we obtain that for 
every t > 

/ (a-'b-H) + a^f [aVt] ^ 2af{t), 

which implies the required inequality. □ 
The proof of Theorem 1 presented below can be essentially shortened by 
using the main result from |T4]. Nevertheless we decided to give an immediate 
and elementary argument following some ideas of [9]. However, first notice 
the following fact concerning subgroups of the multiplicative group (0, oo). 



10 



Remark 1. It is well known (cf. [8], Ch. XXIII, Th. 438) that every 
subgroup of ((0, oo), ■) is either of the form {g": n E Z} with some q G [1, oo), 
or is a dense subset of (0, cxd). The second case occurs, for instance, when 
at least two elements x, y of the group have non-commensurable logarithms: 
In x/ln y ^ Q or the group contains elements different from one but arbitrarily 
close to it. In particular, each subgroup with non-empty interior is equal to 
(0,oo). 

Proof of Theorem 1. By Lemma 3 we may assume, without loss of 
generality, that f{t) > for every t > 0. It follows from Lemma 2 that 
a, 6 G (0, 1). For every t > we define the orbit C{t) of t by 

C{t) = {a'b>t: i,j G Z} . 

Since C(l) is a subgroup of the group (0, oo), Remark 1 yields that either 
In a/In 6 G Q and C{t) = {cH: j G Z} for some c > 1, or In a/In 6 ^ Q and 
the set C(t) is dense in (0, oo). 

Take to > and let G = {lnx:x G C(l)} and A = {lnx:x G C{to)}. 
The function F: A i-^ M, defined by F{x) = ln/(e'^), is continuous and, by 
Lemma 6, convex: 

2F{x) ^F{x-h) + F{x + h), x e A, h e G. 

Therefore, if G = M then, according to [13j VII. 3, Th. 2, the function 

F{x + h)- F{x) 



Ax iG\{0}) 3 {x,h) 



h 



is increasing with respect to each variable. Clearly the same holds true if the 
group G is discrete. In particular, we obtain that the following limits exist 
and do not depend on h: 

Fix + h) - Fix) , , Fix + h) - Fix) 

p_ = lim ; , and p+ = lim 



h x^oo h 

Obviously — oo ^ p_ ^ p+ ^ oo. If p+ = oo, then, taking x = Int, h = \iaa 
and /i = Info, we would have 

lim } = 0, and lim '^^ — ^ 



t— >oo 



fit) ' *-oo fit) 
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which is impossible in view of equation (1). Thus p+ < oo. Similarly p_ > 
— oo. Now we have that 



and 



lim —r--- = aF~ , lim —rrr = oF"^ , 
t->o f{t) t-.^ fit) 



nm^^6^-, Umf^^5^^ 
fit) t^^ fit) 

The equation (1) implies that a^- +6^- = 1 and aP+ +bP+ = 1. As a, 6 G (0, 1) 
there is exactly one real root p of the equation + IP = 1; clearly p > 0. 
Consequently, we get p^ = p+ = p. 

Notice now that the monotonicity of A 9 x ^— > + h) — Fix)) 

implies the monotonicity of the continuous functions C(to) 9 t t— > fi(it)/ fit) 
and C(to) 3 t ^^ and thus 



/(at) = aP/(t), and fibt) = Ffit), t G C(to) \ {0}. 
If t G C(to) then t = a'^b^to for some j G Z, so 

/(t) = /(a^fe^to) = a^'b^^fito) = 

This implies that the positive and continuous function H on (0, oo), defined 
by 

Hit) = t-Pfit), 

is constant on every orbit, and thus, in particular. Hit) = if(at) = Hibt) 
for each t G (0, oo). To complete the proof it is enough to observe that if 
In a/In 6 ^ Q then C(to) = [0, oo) for each to > 0, and thus H is constant on 
[0, oo). □ 

Observe that the property of H in Theorem 1 is similar to that of dou- 
ble periodicity in complex function theory, which goes back to Jacobi and 
Weierstrass (see eg. [4], Ch. XIII, XIV). 

3 Stable distributions in the sense of the weak 
generahzed convolution 



In [25] K. Urbanik considered stable distributions with respect to a gen- 
eralized convolution o. Theorem 4 in [25| states that if /i is a continuous 
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homomorphism from V+ to R and for some measure A £ V+ there exists a 
sequence (c„)ngN of positive numbers such that 

^c„A^" ^ Ao 7^ 5o, 

then there exist positive numbers c and p such that 

h{T^Xo) = exp {-cxP} , x^O. 

The existence of at least one such measure A follows from the property (v) of 
the definition of the generalized convolution o; thus we know that for X = Si 
such a sequence (c„)„gN exists. Because of the form of /(.(T^-Aq) it was natural 
to call the measure Aq stable with respect to the convolution o. 

For the weak generalized convolution the condition (v) need not hold, 
so we will define ©^-stable distributions using the classical linearity condi- 
tions. In this paper we restrict our attention to the symmetric weakly stable 
distributions. 

Definition 2 Let fi be a non-trivial symmetric weakly stable m.easure on 
a separable Banach space E. A measwe A G V+ is stable with respect to the 
weak generalized convolution ® = ®^ if 

V r, s ^ 3 c(r, s) ^ 0, d{r, s) e K (T.A) © (T.A) = (r,(,,,)A) © 5dir,s). 

If for every r, s > we have d{r,s) = then we say that A is strictly stable 
with respect to ©. 

Remark 2. The condition in Definition 2 written in the language of 
random variables states that there are functions c: [0, oo)^ —>■ [0, oo) and 
d: [0, ooY R such that 

V r, s ^ rXe + sX.'e' = c(r, s)X^ + d{r, s)X', 

where £(X) = L{X.') = /x, £ {6) = £{6') = A, X,X',9,e' are independent. 
Putting II (r, s)||2 = y/r^ + s'^ for every r, s ^ 0, (r, s) ^ (0, 0), we have 
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whence 

V r, s ^ rX9 + sX.'9' = ci(r, s)X9 + di{r, s)X', 
where ci and di are given by ci(0, 0) = di{0, 0) = 

c,(r..) = |l(r..)lhc(^,^), 

and 

In the other words, in Definition 2 we can always assume that the functions 
c and d are homogeneous: 

c{ru, su) — uc{r, s), and d{ru, su) — ud{r, s), 

for every r, s,u ^ 0. Notice also that, since X has a symmetric distribution, 
aX = — aX for every a G M. This implies that without loss of generality 
we can assume that d is nonnegative, taking if necessary |d(r, s)| instead of 
d{r, s) for every r, s ^ 0. In what follows we will always do so. On the other 
hand, as we will see in Remark 4, the functions c and d, at least in the case 
of symmetric stable distribution /j, = jp, need be neither homogeneous, nor 
unique, nor continuous. 

Notice that for every a; G R the measure 5x is stable with respect to 
the weak generalized convolution © = for each weakly stable measure 
fi. Indeed, for every r, s > we have {Tr6x) © (TsSx) = (TrSx) © Ssx, so the 
condition holds with c(r, s) = r and d{r, s) = sx. Consequently, we will say 
that 5x is the trivial example of ©^-stable distribution, just as in the classical 
case. 

In order to characterize stable distributions with respect to the weak gen- 
eralized convolution ©^ we choose first ^ G E* such that the random variable 
< ^,X >, jC(X) = yu, is symmetric and non-trivial. Since the distribution 

= C{< ^, X >) is weakly stable for each weakly stable measure we have 
that if /i^ o Xi — fj,^ o A2 for some Ai, A2 G V+ then Ai = A2. 

Following the Urbanik construction for every ^ G E* we define a homo- 
morphism h^: V+ — > M by the formula 

h{y)= I HtiHdt), vev+. 
J [0,00) 
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Let ijj{t) — Ji{t$,), t e R. We see that ijj is the characteristic function of 
the measure fi^ and the function W 3 t h(^{Ttv) is the characteristic func- 
tion of the measure ji^o u. From the previous considerations it follows that 
the functional h^(Tf), t G M, separate points in V+ in the sense that if 
h^{Tt\i) = h^{Tt\2) for some Ai,A2 G 7^+ and for each ^ e E* and t e R, 
then Ai = A2. 

Now for every t G M let ip{t) = h^(TtX), where A G P is stable with respect 
to the generalized convolution By Definition 2 there exist nonnegative 
functions c and d on [0, 00)'^ such that the functional equation 

^{rtMst) = (^(c(r, s)t)ip{d{r, s)t) (2) 

is satisfied. According to Remark 2 we can also assume, if necessary, that c 
and d are homogenous. We need to solve equation (2) in the set $2 of pairs 
of characteristic functions, given by 

$2 = i ((^, V') G $ X $: 3 A G P+ V t G M ^{t) = / ilj{ts)X{ds) 

In what follows we consider four possible, complementary but not disjoint, 
cases concerning the function d: 

3.1 measures strictly stable with respect to the weak generalized convolution, 
when d{r, s) = for every r, s > 0; 

3.2 p-self-decomposable measures with respect to the weak generalized con- 
volution, when such that d{r, 0) > or d{0,r) > for some r > 0; 

3.3 semi-stable measures with respect to the weak generalized convolution, 
when d{r, s) = for some r, s > 0; 

3.4 (c, (i)-pseudostable measures with respect to the weak generalized convo- 
lution, when d{r, s) > for each r,s > 0. 

3.1 Measures strictly stable with respect to the weak 
generalized convolution 

Definition 2 states that if d{r, s) = 0, for every r, s ^ 0, then the measure A 
is strictly stable with respect to (B^. In this case we do not need to assume 
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that the measure /i is symmetric, so (2) leads to the equation 



V r, s ^ 3 c(r, s) ^ Vt G M ip{rt)ip{st) = <^{c{r, s)t). 



(3) 



Equation (3) is the classical one defining strictly stable distributions and 
neither symmetry of the measure n, nor homogeneity of the function c is 
required to obtain that ip is the characteristic function of a strictly stable 
distribution. This means that there exist p G (0,2], a > and /3 G [—1, 1] 
such that 



Since (p is the characteristic function of the random variable < ^,X6'>, 
where C(X.) = /x, C{9) = A , X and 9 independent, then we obtain that every 
one-dimensional projection of the random vector X.9 is strictly stable. This 
implies (for details see e.g. Th. 2.1.5 in [2T]) that the index of stability does 
not depend on ^ and the random vector X.9 is strictly stable. In this way we 
proved the following: 

Theorem 2 Let ^ he a non-trivial weakly stable distribution on E. // 
\ E V is strictly stable with respect to the weak generalized convolution 
then n o \ is strictly stable, i.e. jj, and A are factors of a strictly stable 
distribution. 

Example 1. Let U" be the random vector with the uniform distribution 
11 = Un on the unit sphere Sn-i C M" (as well we can consider here any 
projection U"''^ of U" into M''', k < n). The c^n-weakly Gaussian random 
variable r„ is defined by the following equation: 



where U" and r„ are independent, X is an n-dimensional Gaussian random 
vector with independent identically distributed coordinates. It is known, 
and it was already known to Schoenberg in 1938 (see [23]), that for every 

spherically invariant random vector Y in we have Y = U" ■ ||Y||2, where 
U" and II Y II 2 are independent. This implies that r„ is the distribution of 
II X II 2. Simple calculations show that has the density given by 




■n 



■ Tn — (-^1, • • • , Xn) — X, 



/2,n(r) 



2 



2"/2r(f) 
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For n — 2 this is the Rayleigh distribution with parameter A = 2, thus the 
Rayleigh distribution is cj2-weakly Gaussian. For n = 3 this is the Maxwell 
distribution with parameter A = 2, thus the Maxwell distribution is 
weakly Gaussian. Recall that the generalized Gamma distribution with pa- 
rameters A,p, a > (notation r(A,p, a)) has density function given by 



Thus we have that the generalized Gamma distribution r(A, n, 2) is 
weakly Gaussian. 

Now let 0^ be an a;„- weakly strictly a-stable random variable. Then 
is rotationally invariant a-stable random vector for the vector U" indepen- 
dent of 6^. On the other hand every rotationally invariant ct-stable random 
vector has the same distribution as Yy^5^, where Y is a rotationally in- 
variant Gaussian random vector independent of the nonnegative variable 9a/2 
with the Laplace transform e~*°^^ . Finally we have 



for U", r„ and 9a/2 independent. This implies that the density of a weakly 
strictly a-stable random variable 0'^, is given by 




x > 0. 





In particular, if we take a 



1 then 



/i/2(s) 



1 



-3/2g-l/(2x)^ 



V27r 



X > 0. 



Simple calculations and the duplication formula 




show that 



/l,n(r-) 



22-«r(n) 



r > 0, 



r(n/2)r(n/2) (r2 + l)(n+i)/2 ' 
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is the density function of the cUn-weakly strictly Cauchy distribution. 

Remark 3. Finally we get the following interesting statement: 

If X is a nonnegative random variable with the density function fa,n, or 
if X is a random variable such that £(|X|) has the density fa,n and X' has 
also this property, then the random variable 

Y = ||XU" + X'U'"||2 
has the density 2-i/"/„,„ (x • 2^1/°), i.e. Y = 2V"|X|. 

Remark 4. Spherically invariant (or spherically generated, or rotationally 
invariant) measures, mentioned in Example 1, are extensively studied and ap- 
plied in stochastic modeling. More information about such measures can be 
found in [6|. It is worth mentioning here that in 1963 J.F.C. Kingman (see 
|10| ) constructed an independent increments, two-dimensional (in the sim- 
plest case) stochastic process {Xt:t ^0}, where increments are spherically 
invariant. This process is associated to [Yt = \\Xt\\2-t ^ 0}, describing the 
distance of the particle to the origin, has increments independent in the sense 
of the a;2-weak generalized convolution, but obviously, these increments are 
not independent in the usual sense. This paper was an important part of the 
original motivation for Urbanik's generalized convolution. 

The next theorem gives the full characterization of 7o-weakly strictly 
stable distribution A for 7„ strictly a-stable. Recall that among stable dis- 
tributions only strictly stable distributions are weakly stable and, except the 
symmetric case, only M+-weak stability can be considered, i.e. constants 
a, b, c in definition of weakly stable distribution shall be positive (for details 
see [I7]). Non-symmetric stable distributions are not weakly stable on R. 

Theorem 3 Let 7q, be a strictly a-stable distribution on R". Then the 
following conditions are equivalent: 

1. \ is '-fa-weakly strictly stable; 

2. A = €{06^") for some a > 0, p G (0, 1] and Op is a positive random 
variable with Laplace transform Ee~*^p = e~*^ if p G (0, 1), and 61 = 1; 

3. 7a o A = Ta'jap for some ap-strictly stable distribution ■jap and some 
numbers a > and p G (0, 1) . 
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Proof. By Theorem 2 we have that condition 1 implies that the scale 
mixture 7q, o A is strictly stable; thus 1 implies 3. Evidently 2 yields 1. The 
implication 3^2 follows from [3], where M. Borowiecka-Olszewska proved 
that the scale mixture of strictly a-stable distribution is stable if and only if 
the mixing distribution A equals C{a6l^") for some a > and p e (0, 1]. □ 

3.2 c-selfdecomposable measures with respect to the 
weak generahzed convolution 

The case which we are considering here seems to be rather unrealistic but it 
leads to a very interesting class of distributions and, because of this, is worth 
including. 

Assume, for instance, that for some r > we have d{r, 0) 7^ 0, and so, 
by the homogeneity condition d{l, 0) 7^ 0. If now c(l, 0) = 0, then we would 
have ip(t) = ip{d{l,0)t), t G M, and the corresponding measure A would be 
concentrated at a point. This case is considered as trivial. Thus assume 
also that c(l, 0) 7^ 0. Now, putting a = c(0, 1) and P = d{0, 1), we see that 
a, P > and ip) is a solution of the functional equation 

This recalls the definition of 6-semi-selfdecomposable distributions, which in 
the book of Sato [22] was given in the following way: 

A non-trivial probability measure v is semi-self decomposable if there exist 
a number h> 1 and an infinitely divisible probability measure ph such that 

d{t) = d{b-H)pb{t), t e R. 



Notice that the condition 6 > 1 can be omitted if we are talking about 
nontrivial distributions. If 6 = 1, then the condition holds for any u with the 
trivial measure pi =60. For 6 < 1, notice first that 

u{t) = d{h-H)p,{t) = i)ib-H)ptib-H)pbit) 

n-l 

= ... = i}{b-''t)l[pb{b-H). 

j=0 



19 



Substituting t — > b^t, we obtain 

n-1 

u{b"t)^u{t)l[p,(b"-H). 

j=0 

The left hand side of this formula tends to 1 when n oo; thus there exists 
also the corresponding limit of the right hand side. Since absolute values of 
all the functions here are bounded from above by 1, this equality holds only 
if for every t G M 



= lim 



n— +00 



n-1 

i=o 

which is impossible as u is non-trivial. 

In our case u — (p, b~^ — a and pb = ip{(3-), but we cannot assume that ip 
is an infinitely divisible characteristic function. All we know is that ijj is the 
characteristic function of some nontrivial weakly stable distribution. Notice 
that 

n— 1 00 

i=o j=o 

since we have already shown that a < 1. Now we see that ip is the charac- 
teristic function of the random variable 

00 

j=0 

where Xj, j — 0, 1, . . ., are independent copies of the variable X with dis- 
tribution p. Since the class {/x o A: A G V} of mixtures of the weakly stable 
distribution p forms a set which is weakly closed and closed under convolution 
and rescaling. we see that Y = Xd for some 6 independent of X, C{6) = A, 
under the assumption that the series defining Y converges in distribution. 
Now let us accept 

Definition 3 Let X be a non-trivial weakly stable random vector with the 
distribution p and let a G (0, 1). A random variable 9 (or its distribution) is 
-self decomposable in the sense of tlie generalized weak convolution if 
there exists a random variable Q such that 

= {a9) Q. 
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The previous considerations show that the following proposition holds: 

Proposition 1 Let ^ be a non-trivial weakly stable random vector with 
the distribution /j,, a E (0, 1) and let 9 be a'^ -selfdecomposable in the sense of 
the weak generalized convolution Assume that X and 9 are independent. 
Then 

• if 9 — {a9) (Bix P, then the random variable 9 is uniquely determined by 
the condition 

oo 
j=0 

where Xj 's are independent copies of X. 

• if 9 — {a9) Q for some non-trivial random variable Q, then 9 is 
uniquely determined by the condition 

oc 

X^ = ^q;% Qj, 

j=0 

where Xj 's are independent copies o/X and Qj 's are independent copies 
of Q such that (X.j)j^^ and (Qj) jeN are independent. 

In the case when X has strictly p-stable distribution, we have that the 
series defining Y converges at least in distribution, and 

oo / oo \ 1/p 

Thus the random variable 9 exists and P{^ = (1 — dJ'Y^^^^ = 1, which we 
consider as a trivial solution. 

Let X — Un^i with the distribution u!n,i be the one-dimensional margin 
of the vector U" with the uniform distribution on the unit sphere in M"". 
Then we have Ef/„^i = and VaiUn,! ^ 1 < oo. Thus the series defining Y 
converges in L2, so also almost everywhere. Thus the random variable Y is 
well defined, and there exists 9n such that 

00 

X9n ^ /9 J^a^Xj. 

j=0 
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In the case n — 3 the variable has the uniform distribution on the interval 
[-1, 1], thus 

/ X TT sin (3aH) 

Let G be the cumulative distribution function for Y — Notice that 

for every u G [0, (1 — a)~^] 

G{u) = P {Xi + a {X2 + aX3 + ...)< u/(3} = P {Xi + aY' < u/(3} , 

where Y' is a copy of Y independent of Xi. One can also show that for every 
u e [-/?(! - a)~^, /?(! - a)~^] we have 



Moreover, G as the distribution function of a symmetric random variable has 
the property G{—u) = 1 — G{u) for every m e M. G is also strictly increasing 
on [— — q;)~^,/?(1 — and G{—(3{1 — = 0. Similar conditions 

can be obtained considering C/„^i for an arbitrary n & N. 

3.3 Semi-stable measures with respect to the weak gen- 
erahzed convolution 

Let Z — {(r, s) G (0, 00)^: d{r, s) = 0}. In this subsection we assume that Z 
is nonempty, postponing the case Z — to the next section. We assume also 
that Z ^ (0,00)^, since the opposite case was considered in section 3.1. 

Theorem 4 Let c,d: [0,00)^ [0,00). Assume that Z ^ ^ and Z ^ 
(0,00)^. Let ifji/j-.W — > R be nontrivial characteristic Junctions satisfying 
equation (2). Then there exist p G (0,2] and even continuous functions 
H,K:R\ {0} (0, 00) such that 

= e-l*l^^« and = e"!*™ 

for every i G R \ {0} and 

H{rt) = H{st) = H{c{r, s)t) and K{rt) = K{st) = K{c{r, s)t) 

for every t G M \ {0} and (r, s) G Z . Moreover, 

r^Hirt) + s''H{st) = c(r, sfH^cir, s)t) + d{r, sfK^dir, s)t) 

for every i G R \ {0} and r, s > such that c{r, s) > and (r, s) Z. 
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Proof. Take any (r, s) e Z. Then for every t £ M we have that ip{rt)ip{st) — 
(p{c{r,s)t). The equaUty c{r,s) = would imply that (p{rt)ip{st) = 1 and 
next ifi = 1. Thus c(r, s) > 0, whence ip satisfies the functional equation 

<fi{t)^if{at)ip{bt), (4) 

where a := r/c{r,s) and b := s/c{r,s) are positive numbers. This implies, 
in particular, that ip is the characteristic function of an infinitely divisible 
distribution and, consequently, it does not attain value zero. This property 
can easily be proved directly. First we verify that a.b < 1. Without loss 
of generality we can assume that a ^ b. If a = 1, then by (4), (p{bt) = 1 
for every t > 0, so ip would be identically one. Suppose that a > 1 and 
take any u > such that (p{t) > for every t & {0,u). If t & {0,u), then 
a'H,a~^bt e (0, and thus by (4), 

whence, by induction, 

1 ^ ^p{t) ^ ^{a-n), n e N. 

Now, using the continuity of </?, we deduce that ip{t) — 1 for every t e (0,m), 
and, consequently, for every i > 0. Therefore a < 1 and 6 < 1. Suppose now 
that there exists to > such that (p{tQ) — 0. Then (4) gives (p{tn) — for 
each n e N, where 

tn+i = atn or tn+l = &tn, n G N. 

Since a, 6 < 1 the sequence (tn)neN tends to zero and thus, by the continuity 
of (f, we get (p{0) — 0, which is impossible. 

Now we can define f = — \nip. The function / is non-negative, continuous 
and even, /(O) = 0. Moreover, / as a logarithm of characteristic function 
does not attain value zero in a vicinity of zero. Clearly, / satisfies equation 
(1). By Theorem 1 there exist p G (0, cxd) and an even continuous function 
H:R\ {0} (0, oo) such that for every i e R \ {0} 

f{t) = \t\PH{t) and H{rt) = H{st) = H{c{r, s)t). 

Consequently, 
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Now we show that p and H do not depend on the choice of the point (r, s) G 
Z. Indeed, as H{at) = H{ht) = H{t) for every t e M \ {0}, it follows from 
(4) that dP + V = 1. This and the relation a, 6 G (0, 1) force that p, and then 
consequently iJ, in the representation of / are unique. 

We see now that ip is the characteristic function of a semi-stable dis- 
tribution (see eg. [22], Chapter 3) with the characteristic exponent p, and 
consequently p G (0, 2]. 

Now take any (ro, sq) G (0, oo)^ \ Z . Then we have 

ip{aot)ip{bot) = ip{cot)i^{t), t eR, 

with ao = ro/d{ro,So), bo = So/d{ro,So), Cq = c{ro, So)/d{ro, sq). Since ip is 
positive, so is ip. Define ^f: M ^ M by = — In^/;. The last condition implies 
now that 

g{t) = \t\P (alHiaot) + %H{hot) - c^i/(cot)) , t G M \ {0}, 

so it is enough to define K by K{t) = aQH{aot)+b^H{bot)—c^H{cot). Observe 
that g is positive in a vicinity of zero. Thus, as it satisfies equation (1) like /, 
it follows from Lemma 3 that g is positive. Since g{t) = \t\PK{t), t G M\ {0}, 
also K is positive. To obtain the final assertion it is sufficient to insert the 
forms of (f and into the equation. □ 

Corollary 1 Let c, d: [0, oo)^ — > [0, oo). Assume that d{r, s) > for some 
r, s > and the group generated by the set {r/s: (r, s) G Z} is dense in (0, oo) . 
Let {p,ilj:M. ^ R be nontrivial characteristic functions satisfying equation (2). 
Then there are p G (0, 2] and A,B>Q such that 

ip{t) = exp{-A\t\P}, ipit) = exp{-B\t\P}, t G M, 

and 

A (^rP + s° - c(r, sf^ = Bd{r, sf, r, s ^ 0. 

Proof. Observe that given a function G: M \ {0} M the set 

[r/s: r,s e (0, oo), G{rt) = G{st) for each t G M \ {0}} 

is a subgroup of the group ((0, oo), ■). By Theorem 4 we know that H{rt) = 
H{st) and K{rt) = K{st) for all (r, s) G Z and t G M \ {0}. Thus the 
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density assumption and the continuity of the functions H and K imply that 
H and K are constant and inserting the form of (p and ip into (2) implies the 
assertion. □ 



Remark 3. As follows from Remark 1 the density assumption in Corollary 
1 is satisfied in each of the following cases concerning the set Z of all zeros 
of the function d: 

- interior of Z is non-empty; 

- there are (ri, si), (r2, S2) G Z such that In (ri/si)/ln (r2/s2) is irrational; 

- limn^oorn/ Sn = 1 for some points (r„, Sn) € Z with r„ ^ Sn, n E N. 

Remark 5. Notice that in the case of constant functions H = A and 
K = B we can arbitrarily choose the function c: [0, oo)^ — [0, oo) such that 
c(r, s) ^ ||(r, and define d by 

Remark 6. Recall that a random variable X with the characteristic func- 
tion is semi-stable if there exist constants r,b,c > such that for each 

t e R 

= (ct) e"'\ 

Notice that if X is semi-stable with constants r, 6, c then X + ^(c — 1) is 
semi-stable with constants r, 0, c, thus without lost of generality, we can 
assume that X is semi-stable if this condition holds for 6 = 0. More about 
semi-stable variables, including canonical form of the characteristic function, 
Levy-Khintchine representation and examples, one can find in |22| . 

Sometimes (see eg. [2]) authors assume that r is a natural number greater 
than 1. Such definition do not cover all the possible semi-stable variables, 
but it has a very natural interpretation: there exists r G N and a positive 
constant c > such that 

Xi + ...+Xr = cX, 

where Xi,...Xr are independent copies of X. It was shown in [2] that 
branching processes provides a large class of semi-stable distributions occur- 
ring naturally. 
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3.4 (c, (i)-pseudostable measures with respect to the weak 
generahzed convolution 

Equation (2) was discussed by K. Oleszkiewicz in [H] in the special case 
when /i is a symmetric Gaussian distribution on M and by J. Misiewicz and 
G. Mazurkiewicz (see [16]) in the case when /i is a symmetric p-stable distri- 
bution on R for some p G (0,2). In both cases they have obtained that the 
function ip has to be of the form 

^{t) = <^c,D,g,p(t) = exp {-Cltr - Dm 

for some nonnegative numbers C, D and some g > 0. The problem is that 
not all configurations of these three parameters are available in the sense 
that for g > 2 and some numbers C, D the function (pc,D,q,p is not positive 
definite and, consequently, it cannot be a characteristic function. Otherwise 
we would take Dn for n —>■ oo and then (pc,D„,q,pit) — > exp{— Cltl"^} 
would be positive definite, which, however, is impossible for q > 2. This 
argument shows that for any fixed q > 2 and C = 1 the numbers D have to 
be greater than some positive number. 

Maybe it is more interesting that for some parameters C,D > and some 
q > 2 the function (pc,D,q,p is positive definite. It was shown in [19j and [16] 
that for every r, s ^ 

c(r, s) = II (r, .) II , = (r« + s^f' , rf(r, sY = ^ ( || (r, .) ||^ - || (r, .) . 

Since d is nonnegative we have to require that ||(r, s)||p ^ II ('^j -5)11(7 for all 
r, s > 0, which implies that q ^ p, thus we have the following 

Proposition 2 Let p e (0, 2] and q > 0. If the pair {ipc,D,q,p, 7p) e $2 is 
a solution of equation (2) for some C,D > then q ^ p. 

Nevertheless, for all < p, g ^ 2 and every choice of the parameters 
C,D ^ the function ^c,D,q,p is a characteristic function as a product of a 
symmetric p-stable and a symmetric g-stable characteristic functions. 

In [19] Oleszkiewicz proved that 

• for g G (2, 4] U IJfcL2[^^ ~ 2' ^one of the functions ipc,D,q,2 can be a 
characteristic function; 



26 



• for g G IJfeLi(4^; 4A; + 2) it is possible to find c,d > such that ipc,D,q,2 
is a characteristic function. 

In [IE] Mazurkiewicz and Misiewicz proved that for0<p^l,g>2 one can 
find C,D > such that (pc,D,q,p is a characteristic function. 

Notice that in the situation considered here the information that for fixed 
parameters the function (pc,D,q,p is a characteristic function does not mean 
yet that the pair {(pc,D,q,p, 7p) is a solution of equation (2), since we have also 
the assumption {^c,D,q,p,%) £ '^'2, that is for some A G V+ 

^c,D,qAt) = [ exp {-\tfs^} A(ds), t G M. 
^[0,00) 

Substituting \t\P t and a = q/p we would have 



exp{-Ce -Dt} = 




exp{-tsP} X{ds), t^O, 



which means that the function [0, cxo) 3 t —>■ exp {—Ct"" — Dt} would be 
completely monotonic as the Laplace transform of the random variable 9^, 
where 9 has the 7p- weakly stable distribution A with symmetric p-stable 7^. 

Proposition 3 Let C,D > 0, p e (0,2] and q G U^i((2fc - l)p,2kp). 
Then none of the pairs {^c,D,q,p,%) 'is a solution of equation (2). 

Proof. Let w{t) = Ct" + Dt for every t > 0, where a = q/p. If 
{fc,D,q,p, 7p) is a solution of equation (2) in the class $2 then the formula 
g(t) = exp{—w(t)} defines a completely monotonic function on the positive 
half-line. We see that g'(t) = —w'(t)g(t) < for every t > 0. The second 
derivative takes the form 

g"{t)=g{t) [{w'{t)r-w"{t)]. 

Since a = q/p > 1, w"(t) = Ca{a — for every t > and w'(0+) = d 

then for a G (1, 2) we obtain 

lim g"(t) = —00, 
which means that g is not a completely monotone function. 
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Consider now keN,k> 1, and let a e {{2k-l), 2k). The (2A;)-derivative 
of the function g can be written in the following form: 

where hi{t) = 0, and for every m G N 

m— 1 

hm+i{t) = h'^it) - w'{t) [hm{t) - w^^\t)] , w^^'^t) = C n (a - j) • i"-'"- 

It is easy to see that h2k{t) is a linear combination of elements of the form 
with either Pj ^ a — 2k + 1, or Pj = for every j, thus /i.2fc(0^) is finite 
for a e {{2k - l),2k). We see also that w(2fc)(o+) = +oo, so g^^^\{)+) = 
— oo, which means that g cannot be completely monotonic, contrary to our 
assumptions. □ 



Proposition 4 Let C,D > 0, p e (0, 2] and q > p. If {fc,D,q,p, Ip) e ^2 
is a solution of equation {2), then for each a e (0, 1) the pair {(fic,D,qa,pa, Ipa) 
is also a solution of equation (2) and belongs to $2- 

Proof. It is easy to check that the pair {(pc,D,qa,pa,lpa) satisfies equation 
(2). Thus we shall only prove that the function (pc,D,qa,pa is the mixture of 
the function with respect to some measure from P+. To see this, observe 
first that by our assumptions there exists a measure A e P+ such that 

ex:p{-C\t\'^ - D\t\P} = I exp A(rfs), t G K. 

i[0,oo) 

Let Xp be a symmetric stable random variable with the characteristic func- 
tion R 9 t I— > exp{— Q a random variable with distribution A and 9^ 
a positive random variable with Laplace transform [0, 00) 9 i — > exp{— i"} 
such that Xp, Q and 9a are independent. Then we have 

Eexp {itXX^fQ^/"} = Eexp 
= / exp{-\trsP}X{ds) 

J[0,oo) 

= exp {-C|i|«" - D\t\P^} = ^C,D,qa,pa(t) 
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for every f G R. This shows that ^c,D,qa,pa is a scale mixture of the function 
7^ with respect to the distribution of the random variable Q^^"' , as required. 
□ 

Acknowledgement. The authors are indebted to the referee for his/her 
valuable remarks and comments. 
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